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Abstract 

Electrodynamics of composite materials with fractal geometry is studied in the frame¬ 
work of fractional calculus. This consideration establishes a link between fractal geometry 
of the media and fractional integro-differentiation. The photoconductivity in the vicinity 
of the electrode-electrolyte fractal interface is studied. The methods of fractional calculus 
are employed to obtain an analytical expression for the giant local enhancement of the 
optical electric field inside the fractal composite structure at the condition of the surface 
plasmon excitation. This approach makes it possible to explain experimental data on 
photoconductivity in the nano-electrochemistry. 
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1 Introduction 

It is well known that small-particle composites, like nano-structured noble-metal-dielectric com¬ 
posites, establishes a striking response to the electromagnetic irradiation [H El |3] . This effect 
differs radically from one taking place in both ordinary bulk materials or individual nano¬ 
particles. We are speaking about so-called hot spots, which result from strong localization of 
optical energy on the nanometer scale (less than the optical wavelength). This effect is inherent 
in disordered clusters of nano-particles at the condition of the surface plasmon excitation. 

This effect has a more pronounced manifestation on fractal clusters. Fractals are not transla- 
tionally invariant systems (as any disordered systems), and, therefore, can not transmit running 
waves |1]. As a result, dynamical excitations, like surface plasmons, tend to be localized in frac¬ 
tals. Note, that a superposition of electric helds induced by local surface plasmons results in 
focusing and local giant enhancement of the electric field. This giant enhancement has a ge¬ 
ometrical nature of the scale invariance of fractals, which, at some extend, has the the same 
nature as X-ray diffractive maxima are the geometrical effect of a spatial periodicity of crystals. 

In this paper we employ a concept of fractals and fractional calculus for application to the 
electrolysis enhancement in nano-electrochemistry. An enhancement of the electrolysis current 
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due to a rough metal-electrode interface was first observed in 1926 [5]. The further enhancement 
can be also observed when the electrolysis set up is extend to the light electromagnetic held mm 
[H]. This optical induced current is due interplay between the light and fractional nanostructure 
which is the double electric layer on the fractal electrode interface. This leads to a local giant 
enhancement of electric held. Technically, these “hot spots” lead to the breakdown of the 
double layer capacitor that, eventually, changes the electrode potential (potential of the Stern 
layer). The hot spot phenomenon is well known in optics, and relates to optical properties of 
metal nanoparticles, and is largely due to their ability to produce giant and highly localized 
electromagnetic helds. Two phenomena, closely related to each other, are responsible for giant 
enhancement. The hrst one is localized surface plasmon resonance (see, for examples the reviews 
mm, or monograph m), which is charge density oscillations conhned to conducting (metallic) 
nanostructures. Strong enhancement of the electromagnetic held is possible if its frequency 
is near the frequency of these oscillations m El 0 oroi E]. Another giant enhancement due 
to fractal geometry is a so-called geometrical enhancement [12], obtained analytically in the 
framework of fractional calculus, where a question of its experimental manifestation is still open 

ra . 

Application of fractional integro-differentiation and its link to fractals and fractal geometry 
may shed a new light on studies of properties of composite materials with underlying fractal 
structures [H Eg ESI ETj. The fractal concept makes it possible to involve the application of 
the powerful methods of fractional integro-differentiation, see e.g., [T8lE9|- Extended reviews 
on fractional calculus can be found in [181 EH ED] , and a brief survey is presented in Appendix 
A. 

In this paper, we suggest an analytical consideration, based on fractional calculus, for the 
purpose of obtaining an analytical expression for the electric held in fractal composite media, 
namely in fractal electrode-electrolyte interface. This approach establishes a relation between 
the fractal geometry of the medium and fractional integro-diherentiation, and it is also based 
on our previous studies [T2lE2l[2l]. We suggest a coarse graining procedure for the electric held 
in the Maxwell equation to treat the electric induction term, which is a discontinuous function 
in the fractional composite media. This smoothing procedure makes it possible to involve 
the fractional calculus in (quasi-)continuous media [23] and, in the present particular case, 
to obtain an equation for the electric held in a closed form including the fractal polarization 
charge distribution. The main idea is due a seminal result [23], where a link between fractal 
geometry and fractional integro-diherentiation is constituted in the procedure of averaging an 
extensive physical value. It should be admitted that this description of properties of fractal 
media in the framework of the fractional calculus is possible only in some approximation. This 
procedure of hltering, expressed by means of a smooth function over a Cantor set, leads to 
fractional integration. In its eventual form, it has been presented in Ref. [25]. This topics 
was extensively studied [231 E3] in application for signal processing and dielectric relaxation in 
complex media [2S1 EZl EB] • 

The main idea of hltering, or embedding a matter inside a fractal is the construction of a 
convolution integral according integration over the fractal volume fi{x) ~ x°‘ 

{F(i)> = j F(x)d^(x) ^ „I‘F(x) , (1) 
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where qI^ designates the fractional integral of the order of a in the range from 0 to x: 


„I;F(x) = /V - xr-^F{x')dx'. (2) 

f (tt) Jo 

Here r(a) is a gamma fnnction and a is the fractal dimension. For example, in Ref. [25] 
this kind of integration has been obtained by means of averaging over one of a parameter, 
which describes a fractal strnctnre. Considering a convolntion integral F{x) = W{x) f{x) = 

W{x — y)f{y)dy with the fnnction hF(x), which obeys the scaling relation W(x) = ^W{bx) 
and describes a Cantor set, one performs averaging over the parameter b (namely Info [25]) 
and obtains Eq. ([T]), where the averaging reads {F{x))irib = A{a,b)oI^f{x), where A{a,b) is a 
constant factor and a = |^. Instead parameters a and b one introdnces the fractal dimension 
a and the log-periodicity parameter In 6. Therefore, averaging over the log-periodicity means 
averaging over the all possible realizations of the fractal. 

It shonld be stressed that this idea was hrstly expressed by Kolmogorov in [29] for fractional 
Brownian motion, and later rediscovered by Mandelbrot and Van Ness, who explored this 
phenomenon in greater detail in [30] . 

As shown [121 Ell |22] , these mathematical constrnctions are also relevant to the stndy of 
electrostatics of real composite structures. The main objective of the present research is to 
solve a standard electrostatic problem, namely, the derivation of the electric held in a fractal 
metal-dielectric composite on the electrode-electrolyte interface in the presence of the double 
layer capacitor. This leads to the exploration of the fractional Maxwell equations. This issue 
attracts much attention and is also extensively studied to explore the role of fractional calculus 
in electrodynamics [SU |32l[33l[3ll ES] (see review [36]). 

This formulation can be also employed for the fractal geometry, in particular, for the explo¬ 
ration of the electric held in fractal dielectric composites [37] El], or metal-dielectric composites 
[El El]. This problem is important for nanoplasmonics, where the interplay between the light 
and fractional nanostructures leads to a local giant enhancement of the electric held [HE]. 

2 Electro-mechanical consideration of electron motion in 
the double layer 

The double layer plays a central role for stability in electrolysis current phenomenon near the 
electrodes, see, for example, [SB]- An electron motion inside the double layer in the presence of 
the external high-frequency electric held is described by equation 

mx = —^ -I- f{ujt). (3) 

ax 

Here U{x) is the double layer constant held potential, which is considered as a capacitor po¬ 
tential, while f{ujt) = focosut = eEoit) is a fast-oscillating perturbation, such that uT S> 1, 
where T is an averaged time propagation through the double layer. Here Eo{t) is the electric 
component of the external light and e the electron charge. It is well known that the average 
procedure over the fast oscillations leads to the additional term to the potential due to the 
kinetic term [39] 

= + = + w 
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This additional term leads, eventually, to effective decrease of the double layer capacitor 

potential W —)■ We fcr, where k^T is the Boltzman temperature. To show this let us consider 
the Poisson-Boltzmann theory of the diffusive double layer |38] . 

2.1 The Poisson-Boltzmann equation 

Now we follow the Poisson-Boltzmann theory near the double layer [38]. For the considered 
planar surface, the x coordinate determines the direction of the potential change in the double 
layer. Therefore, the electric potential 'ijj{x) near a charged planar electrode interface with 
the surface potential 'il^i^x = 0) = i/'o, taken as the boundary condition for the double layer, is 
determined by the one-dimensional Poisson equation 

d^ij{x) p{x) 

, 2 =-’ 

dx^ Ed 

where p{x) is the local electric charge density and Ed is the permittivity of dielectric solution 
(at this point it does not carry any important information). The charge density is determined 
by the Boltzmann statistics that reads for the local ion density n* = where Ai is 

the performed work to bring an ion from inhnity to the point x inside the double layer. In 
the presence of the additional potential term t// in (jl|) due to breakdown of the double layer 
capacitor, the electric work Ai depends on the ion charge. Namely, to bring cation the required 
work is A_ — Uf = —e'ip{x), while for anion it is Aj^ + Uj = e'lp- Therefore, the effective potential 
decreases, and since Uf is independent of the coordinates x, this leads to only to re-scaling of 
the boundary condition psi^ —)■ "00 ~ Uf/e. This leads to change of the constant value of the 
solution of the Poisson-Boltzman equation ([5|) 'ijj{x) = {'ipo — [///e)T(a;), where \h(a; = 0) = 1. 

To see this effect of the external light influence on the electrolysis, one needs t// ~ ksT. If 
the electric component of the external light, with u ~ 10^®, is ~ 1 10 V/cm, the additional 

potential energy is 17/ ~ (10“® which is vanishingly smaller than temperature 

fluctuations. Therefore, a condition of observation this effect of the influence of the external 
light on the electrode potential is an essential enhancement of the electric held of the order 
of 10^ 10® times due to the surface plasmon resonance that yields the respond/enhanced 

electric held of the order of 10® 10^(V/cm). Ehect of this giant enhancement of the external 

low-intensity light is considered in the quasi-static approximation. 


3 Quasi-static Electromagnetic Field 

The fractal nanostructure placed on the electrode surface is embedded in the 3D space, there¬ 
fore, we consider a sample of a size R which consists of a fractal metal nanosystem embedded 
in a dielectric (electrolyte) host medium. The system is subject to an external electromagnetic 
held Eo(t) at an optical frequency u, and in the experimental setup, the double layer fractal- 
host inhomogeneities (of the order of R) are much smaller than the light wavelength A. In this 
case, a quasi-static approximation is valid jlO]. Therefore, the nanosystem contains a metallic 
fractal with a fractal volume Vn ~ R^ and the permittivity Em{oo), which depends on optical 
frequency u, while the dielectric-electrolyte host has the volume Vh and the permittivity Ed- At 
the considered scale, the magnetic component is not important and cannot be considered. The 
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electric component satisfies the Maxwell equation in the Fourier frequency domain 


V- 


e(r,a;)E(r,cu) 


0 . 


(6) 


The permittivity can be expressed by means of a characteristic function x(r) [TTl Hll |T2] in 
the form £(r,a;) = +^d[l ~x(r)], where the characteristic function inside the fractal 

is x(r) = 1, r{x,y,z) G Vd, while inside the dielectric host it reads x(r) = 0, r{x,y,z) G 14. 
Following Ref. [12], one splits the electric held into the two components E(r) = E(r) + Ei(r) 
and looks for the reaction of the nanosystem on the external held, where E(r) is the electric 
held induced by Eq in the homogeneous nanosystem, when y;(r) = 0 for V r, or y;(r) = 1 for 
V r and VE(r) = 0. Here Ei(r) results from inhomogeneity of the nanostructute. Thus Eq. 
(j6i) can be rewritten in the form 

x(r)V ■ Ei(r) - g(a;)V ■ Ei(r) + Ei(r) ■ Vx(r) = -E(r) ■ Vx(r), (7) 

where q{oj) = In general case, mixed (Dirichlet, or Neumann) boundary conditions 

are imposed mm- It follows from the boundary conditions that V ■ E(r, cj) = 0 with 
E(r G S', cj) = Eo(ci;), where S means boundaries. This also supposes x(r G S') = 0. 


3.1 Filtering inside a random fractal 

Since the characteristic function x(r) is discontinuous, the next step of our consideration is 
coarse-graining the electric held, which is averaging Maxwell’s Eq. (I7j). This procedure relates 
to integration with the characteristic function J x(r)VEi(r)(il/ and evaluation of the integral of 
the polarization charge density term J Vx(r) ■ E(r)(il/. To that end, let us consider a spherical 
volume of the radius r, such that Iq r < R, where Iq is a minimal size of the self-similarity 
of the fractal volume on the electrode corrugated/fractal surface. In the sequel we will work 
with dimensionless variable t/Iq —)■ r. The electric held does not change at this scaling. A 
fractal mass inside the volume is Ai{r) ~ r^, where 0 < D < 3. Therefore, an average density 
of the metallic phase is of the order of r^~^. For the hltering mass inside the fractal, we 
consider convolution in Eqs. ([T]) and ([2]). We also, reasonably, suppose that the random fractal 
composite is isotropic, he., the (averaged) similarity exponents coincides along all directions. 
This yields that one takes into account only radius r and any changes in the inclination and 
azimuth angle directions can be neglected (this means that the averaged fractal aggregate is 
considered as a sphere with a radially dependent hlling factor r^). Then, we have for the 
divergence 

V ■ Ei(r) = V,.E^,i(r) = V,.Ei(r). (8) 

Using property of the characteristic function, we obtain that the measure of the fractal volume 
Vd is yiVn) = f x(I4))dU = / pcLVd, where dVn = sinOdOdcj) is the elemen¬ 

tary fractional volume in the spherical coordinates and p is a fractal density factor. In what 
follows, we will use /i(r) as the fractal volume by means of the fractal density ~ 4’) 

and for simplicity, violate all constants. Therefore, the fractal volume reads 
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Now, filtering inside the fractal due to characteristic function, which depends only on the radius 
r, yields the integrations 


1 

47r 


x{r)VrE^{r)dV 


x{r')['\/rEi{r')]r'^dr' = ( 5{r' — rj)G'{r')dr' 

rj&Vo 


( 10 ) 


Here we dehne G{r) = r‘^Ei{r) and G'{r) = -^G. 
Following Ref. [12] (Theorem 3.1), we obtain 



-^- [ (r 

T{D - 2) io ^ 


r'f-^G'{r')dr' = oI^-^G'{r). 


( 11 ) 


Therefore, we consider the integration in Eq. fllOl) as the convolution integral with the averaged 
fractal density (r — 


Polarization term Now we estimate the integral E{r')'Vr'X{E)r'‘^dr'. The fractal dust 
Vo at the A^th step of the construction consists of balls with the radius A^v. For example, 
Atv ~ ^ 0 - In the limiting case one obtains Vo = limAr^oo U SS]- The characteristic function 
for every ball is ^(Atv) = 0(r — rj) — 0(r — Vj — An). Differentiation of the characteristic 
function on the intervals [rj, Vj + Aat] yields VrX(AAr) = 5{r — rj) — 5{r — rj — An)■ Therefore, 
for any interval An and at r = r^-, integration with the electric held yields 

rllj^Ni,E{r)r‘^VrX{r) = E(r)r^ - E{r + A7v)(r + An)"^ . 

This expression is not zero in the limit An —)■ 0. Let E{rj) is the electric held outside the ball 
Bn and E{rj + An) denotes the internal electric held. The relation between them, due to Eq. 
(8.2) in Ref. [10] for polarization of a dielectric ball, is 


E{rj + An) = 


3£d 


-E{r,) 


(a^) T 2^^ 

Therefore, the shift for the electric held in the limit Aat —)■ 0 is 

em(uj) - Ed 


E{rj) - E{rj + 0) = E{rj)- 


^m(a^) + 2£rf 

Finally, integration of the polarization charge term yields 


( 12 ) 


E'(r')Vr'X(r')r'^(ir' = ollE{r)r^'VrX{'^) = 


^m(a^) + 


/ E{ry^6{r'-rj)dr'. (13) 


r,GVD 


Again, one obtains the integration of the electric held with the fractal density /r(r') = Xlr eVo 
Vj). This corresponds to the fractal volume ([9]), and, hence, we consider the integration in Eq. 

as the convolution integral of Eq. (ITTll with the averaged fractal density (r — r')^~^. Note 
that this procedure is exactly corresponds to Nigmatulin’s arguments on hltering according Eq. 
(ITi) . where the convolution integral exists from the beginning, while in our case, the convolution 
integral appears owing to theorem of Ref. [12] for integration with the fractal volume/measure 
p(r). Finally, commenting this procedure of hltering inside the fractal, it should be stressed 
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that integration with the fractal characteristic function in Eq. (1T0|) leads to the integration 
with the fractal measure f x(r)G'(r)dr = f G'{r)dfi{r) that results from the physical proper¬ 
ties of x(r). On the contrary, integration with the differential dx{r) in Eq. fll3p leads again 
to the non-zero integration with the fractal measure J x'{r)G{r)dr = f G{r)dfi{r) due to the 
properties of the electric held, namely boundary properties of the electric held on the rough 
metal-electrode-electrolyte interface in the double layer capacitor. 

Eventually, we obtain for the polarization charge term in Eq. ([7|) 


E{r')Vr'X{ry^dr' 


f\r 

r(D -2)J,X 


ry-^E{ry^dr' = + E^]. ( 14 ) 


Here Eq = E ■ f is a projection of the external electric held on the radial direction inside the 
chosen spherical volume of the radius r and 




( 15 ) 


Taking all these arguments and results in Eqs. (JUi and (HI, one presents Eq. ([7]) in the 
coarse-graining form 


oiyy'ir) - q{u)yG'{r)-p{u)oiyyr) = 


2p{uj)Eo 
T{D + iy 




( 16 ) 


The Laplace transform can be applied to Eq. flTbl) . Since G{r = 0) = G'(r = 0) = 0 (the 
electric held of fractal charge density diverges slowly than ^ [m m ED), one obtains for 
(^(s) = £[G(r)] due to Eq. fl26|) (see Appendix): 


= -(17, 

EdS-^ s — q{u;)s^ ^ — p{u;) 

Note that the second term in Eq. (|T 6 ll is q{u)G{r). 

Before arriving at the main result, let us consider the limiting cases. For one 

obtains q{u;) = oo and p{u) = 0. This yields Ei{r) = 0, and the solution for the electric held is 
E{r) = Eq = Er- Another limiting case is |£m(i^)| —t oo. In this case q{u) = 0 and p{u) = —1, 
thus El Eq. Important result here is that permittivity of the mixture is approximately Ed 
that corresponds to the well known result in Ref. [10] [see Eq. (9.7)]. 


4 Surface plasmon resonance 

Now we consider a condition of the strong enhancement of the electric held at the surface 
plasmon resonance, when Re[e(a;)] = —2Ed- This also known as the Frohlich resonance [441145] . 
The permittivity of the metallic nanostructure inclusion at the resonance condition is a complex 
value Emy = -f iE 2 , where £ 2/^1 1 that is described by classical Drude formula (see e.g. 

[401140] j Emy = ^0 — , where Up is a so-called plasma frequency, eq is a high-frequency 

lattice dielectric constant, while the attenuation coefficient 7 is small in comparison with the 
resonant frequency. Therefore, we have for the metallic nanostructure Em{i^) = + i£ 2 , where 

El = Re[£m(i:^)] = eo - and £2 = Im[£„ 7 (c 6 ;)] = 
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At a small detuning from the resonance, when i?e[e:(a;)] = —2£ci + S 2 that corresponds to the 
width of the resonance in the frequency domain, p{uj) reaches the maximal values, that yields 


p(a;) =-1 + —(1-i) ^ —(1-i), (18) 

ze2 

These expressions are inserted in Eq. flTTl) . Before carrying out the inverse Laplace transform 
it is reasonable to simplify the second denominator. We recast Eq. flTTl) in the form 


G(s) = 2p{u)Eo'^ 

k=0 


[q{u)s^ + p{uj)]’‘ 

gfc+4 


( 20 ) 


We take into account that for the scale r ^ 1, the Laplace parameter is small s <C 1, and 
the binomial becomes approximately a monomial. The Laplace inversion of Eq. fl20|) can be 
performed using an expression for the Mittag-LefHer function |16] 




„l-/3 f v-ppsr ^1-/3 f 

1_ / !_ —ds = - _ / e"’' V — ds = V ^ ^ _ 

2m Jc s'^-z 271 Jc ^ s'^’^+P ^ r{uk + P) ’ 

k —0 k —0 


( 21 ) 


where C is a suitable contour of integration, starting and hnishing at —oo and i/,/9 > 0. Com¬ 
paring Eqs. (|2n|l and (|2l|l . (3 = 4 and = 1, one obtains for the electric held 

£’i(r) = 2p{u)EorS(^i^ij\p{uj)r]. (22) 


Since the argument of the Mittag-Leffler function is large, its asymptotic behavior is 

Er 


) ~ n, oc —exp 

p\oj) Ed 


/.. .\ 


(23) 


Eventually, we arrived at the exponential (geometrical) enhancement and giant oscillations of 
the respond electric held due to the fractal geometry of the metal-dielectric composite. 


4.1 Discussion on the geometrical enhancement of the electric field 

As shown, the resulting (enhanced) electric held depends on the parameter p(a;), which is the 
pre-factor and the argument of the Mittag-Leffler function (when r ~ 1) in Eq. ([22]), and 
dehned in Eq. (IT5|l . This parameter is absorbtion efflciency in the electrostatic (quasi-static) 
approximation m- In our case, it describes polarization, and it is obtained in Eq. (1T3|1 under 
evaluation of fractal boundary conditions. One should recognize that the linear (classical |15]) 
enhancement of the electric held is always takes place, due to Eq. (1131) that corresponds to the 
enhancement of the electric held by a small particle and that is rehected by the pre-factor in Eq. 
(120|) . For a fractal small composite of many particles the situation dihers essentially. Here the 
geometrical enhancement of the electric held is due to the focusing ehect of fractal clustering 
rehected by the Mittag-Leffler function and it depends on the argument of the complex value 
of p{uj). 











The exponential (geometric) enhancement takes place when |p(ci;)| 3> 1 and Rep(a;) > 0 as 
in Eq. fl22p . This condition is fulhlled for those frequencies u that are in the vicinity of the 
SPR: Rep(a;) = —2ed + A, where A ~ £2 "C Ed- 

We have to admit that the obtained expressions in the exponential forms are the upper 
bound of the electric held enhancement. The enhanced Ei does not exceed of the order of 
10®(V/cm), otherwise the nonlinear effects become important that violates the linear quasi¬ 
static consideration. This restriction yields r{ed/2e2) < 20, which is a reasonable value for 
experimental realizations. Therefore, we have the light wavelength A ~ 10“^cm, fractal inho¬ 
mogeneity size Iq ~ (10“®cm, and edls2 ~ ojt ~ 5 15, where u is the optical frequency, while 

T is the relaxation time. The latter value determines Iq, as well, which was introduced above as 
a minimum self-similarity size. Note that t < Ts = where is the surface relaxation time, 
while vp is the velocity on the Fermi surface (for free electrons). Therefore, from the condition 
car S> 1, one obtains /q > ^ ~ 10 “^cm. 

5 Conclusion 

It should be admitted that a traditional consideration of photocatalysis is based on the energy 
loss due to plasmon’s decay that leads to generation of hot electrons and holes and suitable 
increase of oxidation-reduction processes. This effect is usually considered under investigation 
of macroscopic, electrochemical manifestation of plasmon resonance |3lll7lll8]. The main reason 
for this follows from the consideration of an individual nanoparticle that yields the maximum 
enhancement of the external held of the order of 10^. Therefore, there is no any measurable 
chemical effects under such high frequency held {u ~ 10^®sec“^) with the amplitude of the 
order of 10^(V/cm), particularly under condition that this held exists only in the vicinity of 
the nanoparticle surface (on the distance of the order of the particle radius). 

Following our consideration in Sec. 2, we were able to show that the manifestation of 
breakdowns of Stern double layer-capacitor due to the hot spots possible at the enhanced 
electric held of the order of 10® 10^(V/cm), and the observed ehect is due to the anomalously 

strong electric held. Nevertheless, it should be stressed that these two diherent ehects can 
simultaneously co-exist in the presence of the fractal self-similarity on the surface of nano- 
structured electrodes. 

Therefore, the main goal of the paper was to develop an analytical method of calculation of 
the giant local enhancement of the optical electric held inside the fractal composite structure. 
It should be admitted that nano-structured noble-metal-dielectric composites are a standard 
object of nanoplasmonics. In our consideration this is a nano-structured electrode - electrolyte 
interface. In this sense, we speak about nanoplasmonic electrochemistry. Owing to the analyt¬ 
ical expressions of enhanced electric held, one can analyze the experimental results in the held 
of photo and electrochemistry [7] . It is worth mentioning that all methods of the calculation of 
the enhancement of the electric held by fractal clusters used so for are diherent approximations 

0 . 

We developed an analytical approach, for description of the wave propagation-localization 
in metal-dielectric nanostructures in the quasi-static limit. The method is based on fractional 
calculus and permits to obtain an analytical expressions for the electric held enhancement. This 
approach establishes a link between fractional geometry of the nanostructure and fractional 
integro-diherentiation. An essential (geometrical) enhancement of the electric held is obtained 
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for the surface plasmon resonance at a certain relation between permittivities of the host and 
fractal metallic nanostructure, when Re£m(i^) = —‘^£<1 with a suitable detuning. 

Important part of the analysis is developing convolution integrals that makes it possible 
to treat the fractal structure. The initial Maxwell equation (|6]) is local, since lo/X 1 and 
space heterogeneity is accounted locally by virtue of the characteristic function. An accurate 
treatment of the fractal boundaries and recasting the Maxwell equation in the form of the 
convolution integrals by accounting fractal properties of the composite, eventually, leads to the 
coarse graining equations, which already take into account the space heterogeneity and nonlocal 
nature of the electric held and polarized dipole charges inside the composite. Therefore, the 
heterogeneity, caused by the fractal geometry, is rehected by the convolution of the averaged 
fractal density and the electric held, according Eqs. ffTT]) and (1TT|) . It is necessary to admit 
that this transform from “local” quasi-electrostatics to the nonlocal, which takes into account 
space dispersion of permittivity e{r) is mathematically justihed and rigorous enough |32]. The 
obtained convolution integrals are averaged values, since the fractal density is the averaged 
characteristics of the fractal structure, and, eventually, it determines the space dispersion of 
the permittivity e{r) of the mixture. 

Summarizing, we have to admit that observation of macroscopic Maxwell’s equations is 
related to averaging of microscopic equations [ID] . This procedure for fractal composite media 
is not well dehned so far, since, according fractal’s dehnition, averaging over any hnite volume 
depends on the size of this volume itself |13]. The main idea to overcome this obstacle is to 
refuse the local properties of equations and obtain nonlocal coarse graining Maxwell’s equations, 
which are already averaged. 
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6 Appendix A: Fractional calculus briefly 

Extended reviews of fractional calculus can be found e.g., in [191 [20l [18] . Fractional integration 
of the order of a is dehned by the operator 

aixfix) = [ f{y){x - yT~^dy , 

r(a) Ja 

where a > 0, x > a. Fractional derivation was developed as a generalization of integer order 
derivatives and is dehned as the inverse operation to the fractional integral. Therefore, the 
fractional derivative is dehned as the inverse operator to qJ", namely aD°f{x) = and 

= aD^^- Its explicit form is 

aD^f{x) = —^ [ f{y){x - yy^-'^dy. 

0!J Jjj 

For arbitrary a > 0 this integral diverges, and, as a result of this, a regularization procedure is 
introduced with two alternative dehnitions of a-D“. For an integer n dehned as u — 1 < a < n. 
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one obtains the Riemann-Liouville fractional derivative of the form 


Jn 

^^D:f{x) = ^air^ix), (24) 

and fractional derivative in the Caputo form 

Jn 

. ( 25 ) 

There is no constraint on the lower limit a. For example, when a = 0, one has q^D^x^ = 
) q;,/ 9 > 0. This fractional derivation with the fixed low limit is also called the left 
fractional derivative. Another important property is D°‘I^ = where other indexes are 

omitted for brevity’s sake. A convolution rule for the Laplace transform for 0 < a < 1 

C[i:f{x)] = s-y{s) (26) 

is commonly used in fractional calculus as well. 
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